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Abstract 

Recently the existence of a random critical line in two dimensional Dirac 
fermions is confirmed. In this paper, we focus on its scaling properties, espe- 
cially in the critical region. We treat Dirac fermions in two dimensions with 
two types of randomness, a random site (RS) model and a random hopping 
(RH) model. The RS model belongs to the usual orthogonal class and all 
states are localized. For the RH model, there is an additional symmetry ex- 
pressed by {Ti, 7} = 0. Therefore, although all non-zero energy states localize, 
the localization length diverges at the zero energy. In the weak localization 
region, the generalized Ohm's law in fractional dimensions, d*(< 2), has been 
observed for the RH model. 
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The study of quantum phase transitions driven by randomness has a long history. In 
1958, Anderson discussed absence of diffusion in random systems [Q. In 1979, scaling 
arguments by Thouless et al. [f2j were further developed |3|]. The dimensionless conductance 
g(L) is treated as the only scaling parameter, where L is the linear system size. When the 
randomness is weak (g 3> 1), the metallic behavior of the d-dimensional system implies 
g{L) = aL d ~ 2 with conductivity a ("the Ohm's law"). On the other hand, when the 
randomness is strong (g <C 1), the wavefunction is exponentially localized, which brings 
g = goexp(—L/^). In the scaling argument, the beta function (3 = d\ogg / dlogL = (3(g) 
plays a central role. In the above two asymptotic limits, the explicit forms are given by 
(3( g ) = (d - 2) + c x /g + 0(1/ g 2 ) (g > 1) and /3(g) = \og(g/g ) (g < 1). It means that all 
states are localized and the quantum phase transition is absent in one and two dimensions. 
However random critical points can appear in two dimensions and much effort has been 
devoted to the study. One of the cases is quantum Hall systems where the time-reversal 
symmetry is broken and they belong to the unitary class In spite of the experimental and 
theoretical implications, since it is beyond the weak-coupling regime, it is still difficult to 
understand the critical phenomena. In this paper, another example of the random criticality 
in two dimensions is studied in detail which has been discovered recently. The properties 
on the critical point ( E = ) has been studied intensively ||. Here the scaling properties 
and the weak localization effect are studied in a critical region ( |-E|->0 ). This is a quantum 
critical phenomena near the random critical line. It brings novel weak localization effect 
('Ohm's law in fractional dimensions'). 

Recently, possible existence of a random critical line in two dimensions was suggested || 
and Dirac fermions with link-type randomness in two dimensions were studied numerically 
The random Dirac fermions were realized on a square lattice by the 7r-flux model 
|| with random hopping. Our model preserves the time-reversal symmetry and belongs 
to the orthogonal class. The zero-energy states do not localize but become critical, which 
can be a prototype of critical states in two dimensions. The density of states is ~ \E\ 
without randomness, and becomes ~ \E\ a in the presence of randomness. The singularity 
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is closely related to the appearance of critical states. Similar phenomenon were found in 



different models ||I0|JTT|| , where the density of states has a singularity and non- localized states 
appear in two dimensions. The exponent a changes with strength of the randomness. It 
implies the existence of the random critical line, which is comparable with other critical 
lines as the Tomonaga-Luttinger liquid in one-dimensional quantum systems. The stability 
of the zero-energy states against the random hopping is due to an additional symmetry 
of the Hamiltonian. The random hopping preserves the symmetry in contrast to the site- 
type randomness. This symmetry appears as a sign change of the Hamiltonian under the 
transformation Cj — ► (—iy x+ ^Cj. The corresponding symmetry in an effective field theory is 
denoted by {H, 7} = with a 4 x 4 matrix 7 One of the possible scenarios proposed in 
H is that all non-zero energy states are localized and the zero-energy states are just on the 
critical point. Unfortunately, since the localization length near the zero energy is very large 
and beyond numerically available system sizes, no direct argument has never been given for 
the scenario. In this paper, we treat Dirac fermions with two types of randomness in the 
light of the scaling. Support for the above scenario is given and weak localization effect near 
the random criticality is discussed. 

We study Dirac fermions with two types of randomness (i) random site (RS) model and 
(ii) random hopping (RH) model. The Hamiltonian is given by 

Hp U re ^ ' ( l)"' ! 'CjtjjCj -|- ^ K C{tijCj -\- ^^CjV^Cj (1) 

i=j±(l,0) »=J'±(0,1) » 

with (i) RS model: t tj = 1, V { = R(Wi) and (ii) RH model: t tj = tji = 1 + R(W 2 ), V t = 0, 
where R(Wys are uniform random numbers between [— W, W\. Although we present data 
with W\ = 2.3 and W2 = 1.0, the qualitative feature does not depend on the strength of 
randomness apart from finite-size effects. 

In the absence of randomness i.e. = 1 and = 0, the model is a tight-binding 
model with half a flux (V flux') per plaquette 0. There are two energy bands on the 
magnetic Brillouin zone [— tt, tt) x [0,7r), which touch at two momenta. Near the two mo- 
menta k* {i = 1,2), where the energy gap closes, the low- lying excitations are described 



by massless Dirac fermions in two dimensions. The effective Hamiltonian is given by 
7~L P ure = 2i / (ix W(yi)[(cr 3 (g)ai)d x + (I ^>a 3 )d y ]^(x), where \&(x) is a four component spinor. 
When the Fermi energy lies at zero energy, that is, all the negative-energy eigenstates are 
filled, the Hall conductivity a xy is ill defined. The sign of mass determines the a xy in the 
continuum theory ]T2|-|rjJ]. There is also the following subtlety in the tight-binding model 



with half a flux (V flux') per plaquette and the next-nearest-neighbor hopping t' [15]. The 
a xy is given by t'/\t'\ and, when t' = 0, the system is on the transition point between states 
with different quantum Hall conductivity. It implies that the zero mode carries non-zero 
Hall conductivity. 

Here we briefly review some properties of the above two random systems. The effective 
field theory of the RS model was discussed by mapping to the nonlinear a model ||16|| . It 



predicts the localization of all states and finite density of states at zero energy. Recently, 
whether the density of states at the zero energy is finite or not for random Dirac fermions, 
is controversial |I7yi5|| . A similar model to the RS model with dilute and strong impurities 
(unitary limit) was also discussed and consistent results with the effective field theory were 
obtained [T^J2^] . All those results suggest that the RS model belongs to the usual orthogonal 
class and standard scaling arguments of the Anderson localization || seem to be valid for 
the RS model. On the other hand, in the case of the RH model, non-localized states were 
discovered at zero energy J7| and the density of states vanishes at the zero energy as ~ \E\ a 
||. It is not only a critical point but also forms a random critical line, since the exponent 
a changes with the strength of randomness. In Ref. 0, the appearance of negative a 
was suggested for sufficiently strong randomness. However, even for strong randomness 
W/t = 1.0, the exponent is still positive and the negative a was not observed ||. Our 
model may be a part of the 'longer' critical line. This is analogous to the massless phase 
of spinless fermions with nearest-neighbor interactions, which is a part of the critical line 
called the Tomonaga-Luttinger liquid. It may be possible to construct models with negative 
a based on our model [^TJ . Divergence of the localization length near the zero energy was 



also suggested numerically. It is not allowed in a standard scenario of the two-dimensional 



Anderson localization. 

Let us first discuss the density of states p(E) = (l/L 2 J2i &{E ~ E{)). We diagonalize the 
Hamiltonian of the RS model for finite squares of size L 2 = 30 2 , 40 2 and 50 2 and ensemble 
average over 10000, 5000 and 3360 realizations is performed respectively. The finite size 
effect is small for the density of states and only the result for L 2 = 50 2 is shown in Fig. 1. 
The result for the RH model is also shown for comparison ||. Finite density of states at 
zero energy is created for the RS model, p(E = 0) ^ 0. It is in contrast to the RH model, 
where the density of states vanishes as ~ \E\ a , p{E = 0) = 0. The difference may be related 
to the presence of random criticality in the RH model, which will be discussed later. 

Next let us discuss scaling properties of the Thouless number g(E) = V(E)/A(E), where 
V(E) is an energy shift obtained by replacing periodic boundary condition with antiperiodic 
boundary condition in one direction and A(E) is a local mean level spacing near the energy 
E. The Thouless number g(E) tells us how the wavefunction is extended in the space. 

Numerical results for the g(E) are shown in Fig. 2. They are shown for L 2 = 50 2 and 
average within an energy window is also performed together with the ensemble average. The 
results suggest that, in both cases, the localization length grows near the zero energy. The 
difference is that the growth is 'singular' for the RH model, which is related to the presence 
of random criticality in the RH model. Although the difference is clear between the RS 
model and the RH model, it is crucial to apply scaling arguments to obtain definite results, 
which will be given below. 

In the following, we assume g = g(L, E, W) = F(y = L/£(E, W)) with the localization 
length £, which means g{L, E, W)'s with different L and E are on a single smooth curve 
g = F(y) using the localization length £(E, W) (scaling hypothesis) [|J f22| . We assume a 
functional form of the £ as 

£{E,W) = \E\K vr >l(\E\,W) (2) 

with a smooth function £(e, W) = l + (,i(W)e + ■ ■ -. The £ n 's are chosen so that g(L, E, W)'s 
with different L and E are on a single smooth curve g = F(y). Here the localization 



length £(E, W) is introduced to define a dimensionless parameter y(= L/£(E,W)) and is 
determined by the scaling hypothesis. It is related to the usual localization length £\ oc (E, W) 
of the exponentially localized wavefunction \ip(x,E,W)\ ~ exp(— |x — x |/^i oc (-E, W)) as 
£,i oc (E, W) = c(W)£(E, W). Fitting our numerical results, we obtain 

= for the RS model (Fig. 3), 

P = -0.75 for the RH model (Fig. 4), 

where W\ = 2.3 for the RS model and W2 = 1.0 for the RH model. This implies that (i) RS 
model: all states are localized, and (ii) RH model: all non-zero energy states are localized 
with the localization length £(E, W) which is diverging as E — > 0, and the zero-energy states 
are just on a critical point. 

Let us discuss the above results in the light of the scaling. Assume that the beta function 
for the Thouless number g obeys the scaling form i.e. (3 = dlogg/dlogL = (3(g). For an 
almost metallic state i.e. g > 1, we expect that A(E) ~ 1/L d and V(E) ~ 1/L 2 due to the 
level repulsion (note that V(E) ~ 1/L for pure systems, since there is no level repulsion). 
Therefore g ~ L d ~ 2 and (3(g) = (d — 2) + c±/g + 0(1/ g 2 ) for large g. In particular, for 
d = 2, (3(g) = ci/g + 0(1/ g 2 ) and g ~ ln(L/^) for large g. We confirmed that the results of 
the RS model can be fitted to this form, which is consistent with usual scaling arguments 
of the Anderson localization |§(see Fig. 3). On the other hand, although states near zero 
energy in the RH model are localized in the thermodynamic limit, they behave as critical 
states due to the large localization length beyond the available system size. Then we can 
expect the g behaves as ~ (L/£) 7 (7 < 0) for the available system size. The results for the 
RH model are consistent with this discussion (see Fig. 4). Note that, since g ~ exp(— L/£) 
for large L/£, there is a correction in L/£ to the above expression which is assumed to be 
g = (L/£) 7 (<?o + 9\L/i + ••■)• We obtained good agreement with this consideration for the 
RH model. 

Here we define an anomalous dimension d* as 

d* = Jim (3(g) + 2 = 7 + 2<2 (3) 
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for states in a critical region in two dimensions. It means g~aL d *~ 2 in the weak localization 
regime (g ^> 1). This is a generalized Ohm's law in fractional dimensions between 1 and 
2. For example, we obtain d*~1.8 in the RH model with W2 = 1.0. The emergence of 
this anomalous dimension d* may be due to the multifractal nature of the zero-energy 
wavefunction in the RH model [|7J. 
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FIGURES 



Fig. 1 Density of states for the RS model (black) and the RH model (white), where 
Wi = 2.3, W 2 = 1.0 and L 2 = 50 2 . A finite width 5 = 0.02 is given to the delta functions, 
although the results do not seriously depend on the small change of 5. The line is guide for 
eyes, which is ~ l^l - 39 . 

Fig. 2 The Thouless number g(E,L) for the RS model (black) and the RH model (white), 
where W x = 2.3, W 2 = 1.0 and L 2 = 50 2 . 

Fig. 3 Scaling function F and the localization length £(E) for the RS model. The 
data are g(E, L) near zero energy with W\ = 2.3, L = 30, 40 and 50 and different symbols 
correspond to different Us. The localization length £(E) is £(E)/\E\P = l + 0.057e — 0.78e 2 — 
0.32e 3 + 0.51e 4 , where f3 = and e = \E\. The scaling function is F(y) = log(y/230)(-0.12- 
0.00060y + 0.0000032y 2 ). The curve beyond the data points is guide for eyes. 

Fig. 4 Scaling function F and the localization length £,(E) for the RH model. The data 
are g(E,L) near the zero energy with W 2 = 1.0, L = 30, 40 and 50 and different symbols 
correspond to different L's. £(E)/\Ef = 1 + 0.47e - 0.0057e 2 + 0.20e 3 - 0.10e 4 , where 
/3 = -0.75 and e = \E\. The scaling function is F(y) = y-°- 15 (0.44 - 0.0097y + 0.000085y 2 ). 
The curve beyond the data points is guide for eyes. 
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